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Quarkonium Physics and a str0 ng from Quarkonia 

Junko Shigemitsu a * 

a Physics Department, The Ohio State University, Columbus, Ohio 43210, USA. 

Recent results in Quarkonia are reviewed, including updates on spectroscopy and a s , and a first look at 
quarkonium annihilation decays. 
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1. Introduction 

Quarkonia serve as good testing ground for 
QCD and for theoretical approaches to QCD, in- 
cluding the lattice, due to their rich spectroscopy 
and decay phenomenology. Furthermore, heavy- 
heavy and heavy- light systems, are playing a cru- 
cial role in probing the electroweak sector of the 
Standard Model. Hence, it is very important that 
we have reliable well tested methods for simulat- 
ing heavy quarks. 

Both four component and two component non- 
relativistic (NRQCD) fermions are being used to 
simulate heavy quarks |IJ§]. Most of the theory 
has been reviewed in previous lattice meetings. 
Only a few comments on recent developments in 
the four component approach will be made here. 
The first formulation of four component heavy 
quarks on the lattice was developed by the Fer- 
milab group To date this approach has been 
implemented at relatively low order. At this level 
one works with the same set of parameters as 
in the clover action for light quark simulations, 
namely the Wilson hopping parameter k and the 
clover coefficient c. For heavy quarks, trunca- 
tion to these two parameters, does not allow the 
coefficient of the p 4 term in a low momentum dis- 
persion expansion, E — M\ + 2 /2M2 + 0(p 4 ) to 
come out correctly. What implications this has 
for M^ adron versus M% adron of composite parti- 
cles, is discussed by Andreas Kronfeld in a poster 
session at this conference M. For hadrons con- 
taining heavy quarks, it is M2 adron that should 
be identified with the physical mass. In order to 
minimize uncertainties associated with M^ adron 
and in fixing k it is important that the p 4 con- 
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tribution in the action be improved. Problems 
with M2 in heavy-heavy and heavy-light hadrons 
using a straight tadpole improved clover action, 
were pointed out a year ago at Melbourne by Sara 
Collins and John Sloan jjj. 

Other four component approaches to heavy 
quarks have now appeared. Alford-Klassen- 
Lepage have a D234(2/3) action that goes beyond 
clover and removes 0(a 2 ) errors at tree level. 
Working with anisotropic lattices, they find good 
relativistic dispersion relation and spectroscopy 
results even for heavy quarks around charm ||. 
The MIT-BU group has implemented a "per- 
fect action" for fermions truncated to couplings 
within a hypercube and also find encouraging re- 
sults for charm Q. 

Due to time constraints, for the rest of the talk 
I will concentrate on quarkonium physics results 
rather than formalism and refer the reader to par- 
allel and poster sessions for more theory. 

2. Brief Overview of Spectroscopy 

The first test of any lattice method for heavy 
fermions, is to sec how well the observed quarko- 
nium spectrum can be reproduced. Figures 1 
& 2 summarize Upsilon spectroscopy results by 
several groups using both quenched (rif = 0) 
and dynamical (n/ = 2) configurations. Tad- 
pole improved clover (Fermilab & SCRI) and non- 
relativistic fermions with improvements through 
0{Mv 4 ) (NRQCD) have been employed. Two 
parameters were fixed by experiment, the bare 
quark mass and the lattice spacing. Fermilab 
and SCRI use the S-P splitting to fix a -1 . The 
NRQCD collaboration determines a -1 from both 
the S-P and the 1S-2S splittings and uses an aver- 
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Figure 2. T Spin Splittings : Symbols have the 
same meaning as in Figure 1. 



Figure 1. T Spectrum 



age in the plots. A first glance at Fig. 1 suggests 
good agreement between simulations and exper- 
iment. A closer look reveals, for instance, that 
for the NRQCD data neither the IP nor the 2S 
state lies on the experimental line. This is par- 
ticularly noticeable for the n/ = results and 
is a reflection of a quenching effect that leads to 
different quantities giving different a~ 1 's for the 
wrong number of dynamical flavors if these quan- 
tities are governed by different characteristic en- 
ergy scales. We will see later, that such details 
are important in precise determinations of a s . 

Quenching effects are also observable in the 
spin splittings of Fig. 2 . They tend to underesti- 
mate spin splittings. Hypcrfinc and fine-struture 
splittings are discussed in more detail in later sec- 
tions. Fig. 3 shows lattice spectrum results for 
Charmonium. 
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3. a s from Quarkonia 

One major goal of lattice gauge theory is to 
contribute towards determination of fundamen- 
tal parameters of the Standard Model. Quarko- 



Figurc 3. Charmonium Spectrum : Symbols 
have the same meaning as in Figure 1. 
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nium studies have led to accurate determinations 
of heavy quark masses and of the strong cou- 
pling constant |]|| . a s estimates are continually 
being refined and improved fl(i(| . We present here 
a status report. 

3.1. Procedure 

There are three steps involved in determining 
a s from the lattice and comparing with other de- 
terminations, (i) Define a suitable a s and mea- 
sure its value, (ii) Set the scale at which a s takes 
on this particular value, (iii) Convert to ajjg. 

For step (i), two popular choices are ap ( P: 
plaquette ) and asF ( SF: Schroedinger Func- 
tional ) |12f| . In quarkonium studies to date one 
has used ap. It is defined through the lxl Wil- 
son loop as, 

4-7T 3 4 

-lnWi A = — ap(— )[1 - (1.19 + 0.07n f )a P ] 

ap agrees through order a 2 with the more famil- 
iar av of Brodsky-Lepage-Mackenzie. It is de- 
fined so that the above relation terminates at a 2 
and there are no higher order corrections. Hence 
the only errors in measuring the value of ap 
come from statistical errors in W\ i which are ex- 
tremely small. In principle one must ascertain 
that nonperturbative, i.e. non power series, con- 
tributions to W% t i are negligible. Several tests 
confirm that any nonperturbative contamination 
would be smaller than one percent ||. 

Most of the work and the source of all the un- 
certainty in a s determinations involve steps (ii) 
and (iii). Quarkonium studies come into play at 
step (ii), when one sets the scale Q = — . by in- 
serting a value for a -1 . Any dimensionful quan- 
tity can be used to determine a -1 . Quarkonium 
level splittings are particularly suitable since sys- 
tematic errors are easiest to control. 
Finite volume : Quarkonia, especially bb, are 
small compared to light hadrons 
Mq dependence : There is no need to carry out 
mass extrapolations like the chiral extrapolation 
in light hadron studies. One works directly at or 
very close to realistic heavy quarks masses. Fur- 
thermore experimentally S-P and 1S-2S splittings 
in quarkonia are very insensitive to Mq in the 
charm to bottom quark region. 



Finite a: Heavy quark actions, such as NRQCD, 
can and have been improved through order a 2 . 
Effects of a 2 errors in the Wilson glue action have 
been estimated perturbatively. Fully improved 
gluonic actions through order a 2 and beyond, are 
starting to be employed as well. 
Quenching or rif dependence : This is the most 
subtle systematic error that needs to be under- 
stood and corrected for. Two questions must 
be settled first; 1. What is n phys \ the relevant 
number of dynamical quark flavors ? 2. How do 
quarkonium splittings depend on m,q Vn ', the dy- 
namical quark mass ? 

Typical gluon momenta px inside an T 
are between 0.5 and l.GeV. The relation 
m u , rrid, m s << pr < m c tells us that n phys ' = 3 
for T splitting physics. This answers the first 
question. It is necessary to ask the second ques- 
tion, because dynamical quarks in current sim- 
ulations are heavier than physical up- or down- 
quarks. Perturbation theory would indicate that 
energy level splittings, AE, depend quadrati- 
cally on rn^ yn ' . However, due to chiral symme- 
try breaking in QCD, AE actually has a linear 
m dyn dependence JljJ. The relevant combination 
of light quark masses is (m u + rrid + rn s ) ■ Since all 
three quark masses are significantly smaller than 
px, one should be able to replace the three un- 
equal mass quarks with three flavors of degenerate 
quarks with m e ff = {m u + + m s )/3 ~ m s /3. 
T physics should not care about the details of 
the m Uy d — m s mass differences. In summary, to 
undo quenching effects, one needs to extrapolate, 
n f -> nf ys - = 3 and m^ n - -> m s /3. 
The fact that one only needs to extrapolate in 
m dyn down to m s /3 and not beyond, highlights 
one of the true advantages of working with the 
T system. Unquenching in systems where one 
is more sensitive to m Ut d — m s mass differences, 
e.g. any system where characteristic gluon mo- 
menta are of order Aqcd or m s will be more 
complicated. Even within quarkonia, the char- 
monium system should have larger uncertainties 
associated with mp ,n than T. 

3.2. Results 

For the quenched, rif — 0, case several groups 
have now extracted a -1 from quarkonium S-P 
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splittings, which can then be used to set the 
scale in ap /_0 ' > [— ]. However, very few rif > 
results exist to date. The most accurate num- 
bers come from the NRQCD collaboration (using 
HEMCGC dynamical configurations) who have 
combined this with their n/ = numbers and 
carried out the extrapolation to n p ^ hys described 
above 1^,0. Recently the SCR1 group has cal- 
culated an 71/ = 2 a^ 1 from the T S-P splitting 
on the same dynamical configurations, however 
using tadpole improved clover 1 have taken 
the liberty to combine their a -1 with quenched 
clover results from the Fermilab group to perform 
the rif —> n phys extrapolation. 

Fig. 4 summarizes results for 1/ap versus 
In(^). This is an updated version of a plot orig- 
inally made by Matthew Wingate |l6| . Only data 
from the last two years are included. There are 
several comments to be made. 

The rif — charmonium data (open boxes) and 
T data (open circles) lie on two separate scaling 
curves. This is due to the fact that the charac- 
teristic energy scale differs between charmonium 
and T. In a quenched simulation, a[Q] runs in- 
correctly between the two characteristic scales. 
When real world data (which incorporates cor- 
rect running of a s ) is used to extract a" 1 , a mis- 
match is found which results in the two scaling 
curves being shifted horizontally with respect to 
each other. A similar phenomenon is found by the 
NRQCD collaboration when they compare a _1 's 
from T S-P and 1S-2S splittings. It is only af- 
ter extrapolating to n phys ' that physics becomes 
independent of which splitting was used ( see be- 
low). 

The rif = 2 configurations used by NRQCD 
and SCRI have am dyn - = 0.01. How close is this 
to the desired m s /3 ? Based on the light stag- 
gered spectrum calculations by HEMCGC |i~7| , 
this m dyn - lies anywhere between ~ m s /2 (from 
the kaon) and ~ m s (from the <ft). Hence there 
are large uncertainties in how far to extrapolate 
in m dyn \ Furthermore, in order to do the ex- 
trapolation, data at another value of m dyn ' are 
required. The NRQCD collaboration has results 
from am dvn - = 0.0250]. The errors in those 
simulations are still large. A 1 to 1.5 a difference 



was found between am dyn - = 0.01 and 0.025. It 
is not clear yet whether true m^ y ™ -dependence 
has been observed. So, the 0.025 data has only 
been used to estimate the error associated with 
m dyn. ^he central value for a p is given by the 
0.01 result and the error is estimated by taking 
the difference between this and an extrapolation 
using the 0.025 data to am dyn - = (0.01)/3. 
The NRQCD results after the rif — > 3 extrapola- 
tion in l/a p are, 



a { " f=3) \8.2GeV] = 



.1948 (29) (11) (37) (S-P) 
.1962 (41) (08) (40) (15-25) 



The three errors correspond to statistical, dis- 
cretization/relativistic and to m dyn ' errors re- 
spectively. The S-P result is the rightmost rif = 3 
(fancy box) data point in Fig. 4 ( the IS - 2S 
number falls right on top of it). The other four 
rif — 3 data points correspond from right to left 
to Fermilab/SCRI, CDHW|| and to two prelim- 
inary NRQCD results from T and charm at lower 
(3, using UKQCD quenched and MILC dynamical 
configurations. 

Finally, step (iii) involves conversion to the MS 
scheme, in order to be able to compare with other 
determination of a s . The conversion formula is, 



(e 5/b Q) 
l + -a£ f) +C 2 (n f )a 



The two loop coefficient C 2 has been calculated 
by Luescher & Weisz for the quenched theory, 
C 2 (rif = 0) = 0.96 §|. In the past C 2 = was 
used. It now makes more sense to use C 2 — 0.96 
for central values even for a^ n f= z \ The differ- 
ence between C 2 = and C 2 = 0.96 will serve as 
an estimate of errors due to uncertainties in the 
a p — Oijjg conversion. 

New results for ot-]js\Mz°\ exist this year from the 
NRQCD collaboration, who have updated their 
previous numbers, and from combining Fermi- 
lab's quenched and SCRFs dynamical T data. 
NRQCD quotes (the Luescher- Weisz C 2 has been 
used), 
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.1175 (11)(13)(19) (S-P) 
.1180 (14) (14) (19) (15 - 25) 
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Figure 4. 1/ap versus ln( — ). rif = 
charm: boxes [18,19,20,22,23]; n f = T: cir- 
cles [18,19,21,14]; rif = 2 charm: + [23]; n f = 2 
T: x [9,14,15]; n/ = 3 (extrapolated): fancy 
boxes. Non lattice determinations of a s typically 
find values for ajjg[M z o] between 0.115 ~ 0.125. 

Converting to cJp f ~ 3 ^ , this corresponds to the re- 
gion between the two full curves in the plot. The 
dotted lines are nf = three loop scaling curves. 



The first error now includes both statistical and 
discretization/relativistic errors, the second is 
due to rriq yn ' and the third comes from conversion 
uncertainties. The combined Fermilab/SCRI S-P 
data give, 

a Ws[ M z°} = 0.H59 (19) (13) (19) 

This is a very preliminary number based on this 
author's extrapolations. Nevertheless it is en- 
couraging that results using two different fermion 
methods are approaching one another. 
In the future, improvements in a s from quarko- 
nia should come from further studies of m^ vn ' 
dependence. Simulations with dynamical Wilson 
configurations or staggered with rif > 2 would 
be welcome. The largest source of uncertainty 
for a-jjg- is still higher orders in the conversion 
formula, so Ci(nf > 0) is crucial. Finally, un- 
quenched lattice a s determinations other than 
from quarkonia are eagerly awaited. 

4. Quarkonium Annihilation Decays 

Annihilation decays of quarkonia provide 
another fertile arena for lattice investigations. 
Work has just begun in this area with first re- 
sults coming from the Argonne group p7| . The 
theory was developed a couple of years ago by 
Bodwin-Braaten-Lepage, using the framework of 
NRQCDQ. Annihilation decay rates T(QQ -> 
X) ( X : light hadrons, 77, l + l~), can be ex- 
pressed as a power series in 1/Mq. The coef- 
ficients factorize into a short distance perturba- 
tive and a non-perturbative part. Lattice meth- 
ods can be applied to obtain the non-perturbative 
contributions that enter as hadronic matrix ele- 
ments of four fermion operators, e.g. 

U x = (^1^(1/2) D X-X f (V2) 6 ^P) 

The subscript 8 in Hs denotes a "color octet" 
contribution, which is sensitive the the QQg com- 
ponent in P-wave quarkonia, where the heavy 
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quark and anti-quark form a color octet combi- 
nation. The remaining three singlet contribu- 
tions can be written in terms of quarkonium wave 
functions at the origin or derivatives acting on 
wave functions, using the vacuum saturation ap- 
proximation. For instance, 0i/|( 1 5o|V' x|0)| 2 = 
(l + 0{v 4 )) leads to Q x ~ £\Rs{0)\ 2 ( R ■ radial 

wave function) and Hi/|( 1 Pi|-0 t: ^ D x|0}| 2 = 
(1+C(u 4 )) to Hi - ^^(O)! 2 etc. The Argonne 
group finds that vacuum saturation is well satis- 
fied. They work with the lowest order ( 0(Mv 2 )) 
NRQCD action in their simulations. Their re- 
sults include one-loop matching between lattice 
and continuum operators. A comparison between 
lattice results and experiment is given in Table 1. 
Within large errors, one sees good agreement for 
charmonium, less so for T. Since these matrix el- 
ements are sensitive to wave functions at the ori- 
gin, one expects significant quenching corrections. 
One should also note the large uncertainty coming 
from perturbative matching. The Argonne group 
is looking into the possibility of nonperturbative 
renormalization of the relevant operators. 

5. Spin Splittings 

Spin splittings in quarkonia are sensitive to ev- 
ery single detail in the action. They are signif- 
icantly affected by quenching and have a strong 
dependence on Mq . Tadpole improvement of the 
action, to remove lattice artifacts, was found to be 
crucial for achieving agreement with experiment. 

In Fig. 5 we show charmonium hyperfine split- 
tings in the quenched approximation versus a 2 . 
There are three data points from Fcrmilab us- 
ing heavy clover at /3=5.7, 5.9 and 6.l|ll|. The 
NRQCD point comes from simulations at p = 
5.7^2|. Then there are three coarse lattice re- 
sults using D234(3/2)[^5| on an anisotropic lat- 
tice plus data from Howard Trottier^6|, who ap- 
plied NRQCD on 0(a 2 ) improved glue. The 
"x"'s differ from the "+"'s, in that the former 
have a 2 errors removed from the E and B fields 
( this is an a 2 improvement of the NRQCD ac- 
tion that goes beyond tree-level, and appears to 
be noticeable for a > 0.2/to). It should be 
noted that current NRQCD simulations have er- 
rors of order M c v\ ~ ihMeV for cc and of or- 
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-2.0 — 1.6 GeV 2 
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0.00262(3)(24)t 5 jGeV~ 2 
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Table 1 

Continuum renormalized decay matrix elements 
from lattice calculations, compared with those 
extracted from experimental decay rates where 
available [27]. Errors in the lattice results corre- 
spond from left to right to statistical, systematic 
and perturbative matching errors. 



der MbV® ~ 5MeV for bb. Very preliminary es- 
timates of some of the Mv 6 corrections in char- 
monium indicated the hyperfine splitting could 
come down significantly. Hence, even though the 
current NRQCD number at p = 5.7 is closest 
to experiment, the true quenched result could be 
much lower, e.g. at about the 70MeV level seen 
by the Fermilab group. It is important that all 
the spin dependent Mv 6 corrections be included 
in a future NRQCD simulation. 

6. Summary 

Quarkonia provide many opportunities for lat- 
tice gauge theory. There is a rich phenomenology 
of spectroscopy, decays and yet to be discovered 
states. T S-P and 1S-2S splittings provide well 
controlled scale determinations that lead to ac- 
curate determinations of a s . Spin splittings still 
pose challenges, but we should be up to them. 
Quarkonia serve as good testing ground for new 
actions. Studies of quarkonia create invaluable 
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Figure 5. Quenched Charmonium Hyper- 
fine Splittings versus a 2 . Fermilab : boxes; 
NRQCD : diamond; D234(2/3)3:l : fancy boxes; 
Trotticr :+ & x ; ; experiment : burst. Errors 
are statistical. 



experience for going onto heavy-light systems. 
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